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ABSTRACT 


Pr ev iou s work on efficiency of light injection into the core of a fiber from a 

!u m r-T, and , a bulk distribution of sources in the cladding, have made use of 
t e fields of a weakly guiding fiber. This approximation simplifies the 
analysis of the power efficiency by introducing universal values for the 
eigenvalues of different fibers with same K-number, but can not predict 
accurately the behavior of the injected light into a fiber with arbitrary 
differences in indices of refraction. We have used the exact field solution in 
the expressions of the power efficiency, P eff , and analyzed its behavior as a 

function of the remaining parameters. Although more complicated and 
harder to interpret, our formulas allow us to analyze the power injection 
efficiency of fibers with arbitrary differences in the indices of refraction The 
results obtained are relevant for the design of more efficient optical fiber 
distributed sensors. The conclusions follow. 

We have confirmed weakly guiding results obtained previously. However 
we have found that the P eff does not always increase with the T-number 
but with the difference in the indices of refraction n -n , 

' core ‘clad* 


For fixed a IX, indices of refraction, n core and « clad and normalized inner 
and outer radius, R in and R oui , the P e ^ is independent of the core radius, a, 

and the wavelength, X, for any uniform cylindrical distribution of cladding 
sources. This suggests that a/X, R in and R out are independent variables. 


For the bulk distribution we have found that P eff increases with the 


wavelength, X, and decreases with the fiber core radius, a, 
with the V-number. However, for the thin film, the P e{( ' 


i.e., it decreases 
remains almost 


constant with the wavelength, X, and fiber core radius, a. 


xi 
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INJECTION EFFICIENCY OF BOUND MODES 


I. INTRODUCTION 


1-1 Justification 

In this work we have obtained expressions for the power injection efficiency, 
Prff/ from cladding sources of a step index profile cylindrical optical fiber 
using the Maxwell's Equations. We have also determined the general 
behavior of P eff against the parameters of this problem extending previous 
work done by Marcuse [1988] on the efficiency of core light injection from 
sources in the cladding. Our results are unique in the sense that we have 
made use of the exact field solution of the cylindrical optical fiber to calculate 
the power efficiency due to the emission of evanescent waves and analyzed 
its behavior as a function of parameters and their combinations not 
previously treated. Previous results that are related to our work, have either 
made use of the geometric optics theory to calculate and optimize the output 
signal of a fluorescent coated fiber [Glass et. al., 1987; Love et. al., 1988] or used 
the weakly guiding approximation of the fields to calculate the power 
efficiency [Marcuse, 1988]. The first approach, geometric optics, is accurate 
only for multimode fibers and the second approach, weakly guiding 
approximation, can not predict the behavior of a fiber that has a large 
difference between the core and cladding indices of refraction. In this work, 
we have used the exact field solution in the expression of the power efficiency 
in order to take into account the shortcomings of both approaches. In this 
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way, we have been able to determine accurately the behavior of a few mode 
optical fibers and the behavior of fibers that have an arbitrary difference of 
indices of refraction. Our results can be used to model optical fiber chemical 
sensors that have a fluorescent or chemiluminescent cladding. 

1-2 Overview 

In this Chapter we present an overview of previous theoretical and 
experimental results that are related to this work. In Chapter II, we derive the 
fields of both arbitrary and weakly guiding fibers with an infinite cladding. A 
more detailed presentation of the results in that Chapter can be found in 
many good optical fiber textbooks [Snyder et. al., 1983; Marcuse, 1974; 
Midwinter, 1979] and earlier papers on this subject [Snitzer, 1961; Gloge, 1971 
and references therein]. In Chapter III, we derive the power injection 
efficiency from sources in the cladding using the exact field solution of an 
optical fiber. For the first time, this derivation allows us to predict the 
behavior of fluorescent cladding fibers with arbitrary differences in indices of 
refraction. The weakly guiding approximation referred to earlier, can be used 
only with fibers that have similar core/cladding indices of refraction. 

Chapter IV is central to this work. It presents the numerical results obtained 
from our FORTRAN program which has made use of the formulas derived 
in Chapters II and III. In that Chapter, we have plotted the behavior of the 
power efficiency of a fiber as a function of its parameters namely, the core 
radius, a, the indices of refraction « core and n C | ad , and the wavelength, A. 
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We have also plotted the power efficiency against a combination of these 
parameters as n core -" clad and the V-number under certain conditions. 
There, the results for a new independent variable, a/X, is discussed. Our 
work has confirmed previous weakly guiding results, given a new insight 
into the behavior of an optical fiber with an arbitrary difference in the indices 
of refraction and helped determine the configuration of a highly efficient 
fiber. These results are relevant in that they point toward the design of 
improved optical chemical distributed sensors. 

Finally, in Chapter V, we have concluded and summarized the results 
obtained here and recommended further additional work that could both 
verify and improve our theory. 

1-3 Theoretical and Experimental Background 

Optical fiber sensors have been used to measure current, pressure, 
temperature, and other physical observables [Culshaw, 1984 and references 
therein]. They can also be used to monitor vibrations in flexible structures 
and detect chemical species [Rogowski et. al., 1988; Hardy et. al., 1975, 
Lieberman et. al., 1988; Lieberman et. al., 1990; Blyler et. al., 1988; Wolfbeis et. 
al., 1988]. Many optical techniques have been employed in conjunction with 
optical fibers. They exploit both the wave and corpuscular nature of light. 
Geometrical optics theory is used whenever the wave aspect of light can be 
disregarded and for this reason it is an incomplete theory. 
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Here, we will be mostly concerned with the wave aspect of light. By using 
Maxwell’s Equations, phenomena like evanescent waves can be studied and 
predicted more accurately. It is the concept of evanescent wave that is central 
to this work. It comes into play whenever radiation is totally internally 
reflected between two dielectric media of different indices of refraction. 
Although most of the incident power is reflected, part of the radiation 
penetrates a very thin layer of the dielectric of lower index of refraction. This 

penetrating radiation is termed the evanescent component of the field and is 
characterized by 

1) an amplitude that has an exponential-like decay from the dielectric 
boundary and 

2) a direction of propagation that is parallel to the dielectrics interface 
[Jackson, 1975; Ghatak et. al., 1978]. 


As consequence of 2, all the power in the dielectric medium of lower index 
propagates parallel to the interface (Figure Tl). 

Evanescent waves are responsible for the Goos-Hanchen Effect. It occurs 
whenever a totally internally reflected light beam in a plane dielectric 
interface undergoes a lateral shift due to the penetration of the field into the 
lower refractive index medium or evanescent region [Ghatak et. al., 1978; 

Midwinter, 1979]. The larger the wavelength of the radiation the larger the 
lateral shift. 
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Figure 1-1. The field amplitude of the evanescent wave decays exponentially 
further from the media interface. Its direction of propagation is parallel to the 

interface. 9 C is the critical angle. 
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Evanescent waves can also give rise to other unusual phenomena. They can 
be absorbed and used to excite fluorescent sources in the evanescent region 
[Tai et. al., 1987; Glass et. al., 1987; Carniglia et. al., 1972; Lee et. al., 1979; Love 
et. al., 1988]. The principle of reciprocity, in conjunction with the above 
results, would lead us to expect that the inverse process should also take 
place [Carniglia et. al., 1972]. In other words, if an excited molecule placed near 
the interface of two dielectrics in the lower index of refraction region can emit 
evanescent waves, we should be able to detect radiation propagating at an 
angle higher than the critical angle (Figure 12). Such high angle propagation 
was observed by Selenyi [1913] and Frohlich [1921] and quantitatively 
explained by Carniglia et. al. [1972]. 

In conclusion, evanescent waves can be used to excite sources in the 
evanescent region. The inverse process, emission of evanescent waves by 
excited molecules can also take place and is a well-established phenomenon. 

1-4 Previous Work on Evanescent Wave Sensor 

The evanescent wave concept has been theoretically and experimentally 
investigated and widely used for sensing purposes [Tai et. al., 1987; Glass et. 
al., 1987; Carniglia et. al., 1972; Lee et. al., 1979; Love et. al., 1988; Selenyi, 1913; 
Frohlich, 1921; Lieberman et. al., 1990; Cox et. al., 1985; Lieberman et. al., 1988; 
Christensen et. al., 1989; Marcuse, 1988]. 
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Figure 1-2. Source placed near the interface of two media in the region of 
lower refractive index. Light rays observed at 9>6 C (arrows in medium 1), 
are due to the evanescent waves emitted by the source in medium 2. 
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Tai et. al. have used solely the absorption of the evanescent wave to 
determine the concentration of methane-gas with a tapered optical fiber [Tai 
et. al., 1987]. In their approach, a He-Ne laser excites bound modes in the fiber. 
The chemical species surrounding the tapered region of the fiber, absorbs the 
evanescent wave associated with these modes at a specific wavelength. This 
absorption can be detected at the end of the fiber as a decrease in the output 
signal level and is related to the concentration. 

The effect of the numerical aperture (N.A.) on signal level of multimode 
fibers with fluorescent substances in the evanescent zone was treated (2). 
Using a N.A. smaller than the maximum N.A. of the sensor, the authors 
have excited fluorescent sources in the evanescent zone. Due to the emission 
of evanescent wave, some of the light is trapped in the core as bound modes 
and can be detected at both ends of the fiber. Due to the background noise, the 
signal was chosen to be collected at the proximal end of the fiber. The theory 
developed was shown to be in good agreement with the experimental results. 

The effect of the numerical aperture (N.A.) of injected light on signal level of 
multimode fibers coated with a fluorescent substance in the evanescent zone 
was treated [Glass et. al., 1987]. Using a N.A. smaller than the maximum 
N.A. of the sensor, the authors injected light into a multimode fiber. The 
evanescent field of the radiation excited the fluorescent sources. Due to the 
emission of evanescent wave, some of the light is trapped in the core as 
bound modes and can be detected at both ends of the fiber. Due to the 
background noise, the signal was chosen to be collected at the proximal end of 
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the fiber [Hirschfeld et. al., 1984; 1984a; Block et. al., 1984; Block et. al., 1986]. 
The theory developed used geometric optics and was shown to be in good 
agreement with the experimental results. 

Using the emission of evanescent photons, Lieberman et. al. have developed 
an optical fiber sensor with a fluorescent cladding to detect molecular oxygen 
[Lieberman et. al., 1990]. In this case, an optical fiber is clad during 
manufacture with a polymer which has a fluorescent dye dissolved in it (the 
dye itself is sensitive to the presence of molecular oxygen [Cox et. al., 1985]). 
The dye dissolved in the polymer. Poly-dimethyl Siloxane, PDMS, acts like a 
fluorescent cladding when excited by side-illumination at a wavelength 
within the excitation range of the dye. In a way similar to the one described 
before, evanescent photons from the cladding are injected as bound modes in 
the core and the homogeneous photons are leaked out. In a similar sensor, an 
oxygen sensitive fluorescent coating was applied to a fiber having a 
fluorescent core. The homogeneous photons from the coating sources excited 
the fluorescent sources in the core [Lieberman et. al., 1988]. The result was a 
100-fold increase in the efficiency of the sensor when compared with the 
previous one. In both cases, the fluorescence intensity is quenched by oxygen 
and its partial pressure could be determined [Cox et. al., 1985]. 

Christensen et. al. have used a FDTD (Finite Difference Time Domain) 
analysis to determine the percentage of power that is guided in a single 
direction in a planar waveguide with fluorescent sources outside the core 
[Christensen et. al., 1989]. His approach differed from Marcuse's who derived 
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an expression for the power efficiency of a*step profile fiber using the weakly 
guiding approximation [Marcuse, 1988]. Marcuse has applied his results to the 
sensor described by Lieberman et. al. [1990]. Agreement between theory and 
experiment was within an order of magnitude for a single experimental data 
point. Further experimental work is required to verify the theory. 

In this thesis we extend previous work done by Marcuse [1988] on the 
efficiency of core light injection from sources in the cladding. 

Assuming that the cladding sources are already excited, Marcuse computed 
the fraction of the power that is trapped in one direction as bound modes in 
the core of a weakly-guiding step-profile fiber. This fraction is the ratio of the 
power that is injected into the core of the fiber as bound modes, P core to the 
total power radiated by the sources, P ra j. He also assumed a fiber with 
infinite cladding and treated it for both positively and negatively guiding 
fibers, « cor ^>« c i a d and ti clad > « core respectively. Assuming each source as 
being an infinitesimal electric current j with random phase and orientation, 
he analysed two different distributions: 

1) sources that are uniformly distributed in the cladding and 

2) sources that are concentrated in the core-cladding boundary (thin 
film) (Figure 1-3). 

In this extension, we treat only positively guiding fibers for the above source 
distributions. We derive our formulas by using the exact solution of the fields 
of the cylindrical optical fiber with an infinite cladding. Such an approach 
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(b) 


Figure 1-3. Sources uniformly distributed in the cladding (bulk distribution) 
(a) and sources distributed at the core/ cladding interface (thin film 
distribution) (b). 
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allows us to account for arbitrary differences between » core and « c i ac j. 
Although the cladding sources also excite low loss leaky modes that may be 
detected in short fibers [Marcuse, 1988], we have treated only the injection due 
to the evanescent field. We are applying these results to design an atomic 
oxygen optical fiber sensor which uses evanescent-wavecoupling. Such a 
sensor would be part of a smart structure and could detect and monitor the 
action of atomic oxygen on space structures. It could easily be deployed on 
NASA's missions. 

1-5 Conclusion 

Since optical fibers were initially developed for communications purposes, 
there had been a need to design fibers with low pulse dispersion. Weakly 
guiding fibers seem to fulfill this need. This has led to the formulation and 
widespread usage of the weakly guiding approximation which is simpler than 
the exact case. In contrast, it was only recently that the potential of optical 
fibers for sensing purposes was realized. In consequence, it was also realized 
that weakly guiding fibers are not necessarily best suited for sensing purposes. 
Although the exact field solution can explain sensor fibers a simpler and 
different approach that can treat these fibers has to be formulated. The 
possibility of formulation of a new, innovative and simpler approach can not 
be discarded. It should allow a fair description of the properties of fiber 
sensors in the same way that weakly guiding approximation is a good 
approximation for communication fibers. 


II. BOUND MODES OF A STEP PROFILE FIBER 


II-l Introduction 

The system to be studied, the cylindrical optical fiber, consists basically of a 
glass-rod of index of refraction « core surrounded by a dielectric material of 
index of refraction « c | ad .where 


^core^ ^clad* 

The region with index of refraction n core is referred as the core of the fiber. 
The outside region with index of refraction n c lad is referred as cladding 
(Figure II I). 

A cylindrical optical waveguide can support two kinds of modes: radiation 
and bound modes. The total electromagnetic fields of the waveguide can be 
expressed as a sum of these two modes. In ray theory, light beams that have 
an incident angle smaller than the critical angle correspond to the radiation 
field. Those are the refracted rays. On the other hand, electromagnetic fields 
that propagate indefinitely inside the waveguide structure are expressed as 
bound modes. In general, most of the bound mode energy propagates inside 
the core of the fiber. The portion that penetrates into the cladding, the region 
of lower index of refraction, is termed the evanescent field. As discussed 
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Figure II-l. Three layers cylindrical optical fiber with core radius a and 
cladding radius b (a). Variation of the index of refraction ti with the radius 
r (b). For practical purposes the cladding radius b is assumed to be infinite. 
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previously, there is no ray counterpart for the evanescent field. However, for 
those fields that are not evanescent their geometrical counterpart are 
represented by rays that are totally internally reflected, having an angle of 
incidence greater than the critical angle (Figure II 2). 

Bound and radiation modes are excited by injecting electromagnetic energy 
from sources into the fiber. Lasers, diodes and fluorescent molecules can be 
used for this purpose.. As discussed in the previous chapter, we will be 
concerned only with the excitation of an optical fiber due to fluorescent 
sources distributed in the cladding. 

Excited fluorescent sources distributed in the cladding of an optical fiber 
generate radiation fields and can inject bound modes (Figure II 3). For sensing 
purposes, bound modes are more important. They propagate indefinitely in 
the core of the fiber and can be easily collected for analysis. In addition to that, 
any perturbation to the trapped field can be used to get information on the 
surroundings of the fiber. The bound modes excited by the sources in the fiber 
cladding are closely related to the evanescent field. Without evanescent 
fields, bound modes can not be excited from sources in the cladding. 

It is the aim of this chapter to present a mathematical expression describing 
the bound modes of a cylindrical optical fiber. We start by presenting a well- 
known relationship between the transverse and longitudinal fields. We 
proceed with the derivation of the vector wave equation and find the 
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Figure II-2. Bound modes (rays) and their evanescent counterpart in an 
optical fiber. 
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radiation modes 
-► bound modes 




Figure II-3. Sources in the cladding of an optical fiber can emit evanescent 
wave. Some of the evanescent energy is trapped as forward and backward 
propagating modes. However, in general, most of the it goes toward radiation 
modes. 
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solution of the longitudinal components t>f the fields. From the previous 
results we derive the transverse components of the fields. Using the 
boundary conditions and an appropriate normalization, we find the 
amplitude coefficients and the eigenvalue equations. We conclude the 
chapter by displaying the fields and the eigenvalue equation of a weakly 
guiding fiber. 

1 1-2 Relation Between Transverse and Longitudinal Fields 


Propagation of light along non-absorbing optical waveguides in regions 
sufficiently far from any source of excitation is most conveniently described 
by a sum over the bound modes of the waveguide. Bound modes are 
solutions of the source-free Maxwell equations. We assume an implicit time 
dependence exp(-icot) in the field vectors. The dielectric constant e{x,y,z) 
can be written as a function of the index of refraction n{x,y,z), 

£ = € 0 n 2 . (II-2) 

For nonmagnetic materials the magnetic permeability n is very nearly equal 
to the free-space value fd Q . Under these conditions and using the rationalized 
MKS units, the source-free Maxwell’s equations are expressible in the form 


V x E = H 

V Co 


(II-3a) 
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V xH = -iJ~^fkn 2 E 
V fAo 


V.(n 2 £) = 0 


V.H = 0 


(H-3b) 

(II-3c) 

(H-3d) 


where k = 2k /X is the free-space circular wave number and X is the 
wavelength of the light m free space. 

The refractive index profile of the step index fiber doesn't vary with the 
distance along the waveguide, i.e., it is independent of z. Such fibers are said 
to be translationally invariant [Snyder et. al., 19831- Its electric and magnetic 
fields can be expressed as superposition of fields with the separable form 

E (x,y,z,) = e (x,y ) exp(t 0z) (II ' 4a) 


H (x,yj) = h (x,y ) exp(i f3z) 


(n-4b) 


where (5 is found from the boundary conditions. 

It is useful to separate the fields in components parallel and transverse to the 
fiber axis 


E ( x,y,z ) = (^t (x,y ) (x,y )) exp(ipz) 


(II-5a) 
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H (x,y,z) = (h t (x,y ) + zh z (x,y ) ) exp(/ (h). (II-5b) 

Substituting equations II-5 in II-3 and separating the resulting equations into 
transverse and longitudinal components, we get 


e,= 'V?r L 2 ix ^' ,, + i 

y to k n 2 

( n-6a) 

^ t= V + 

(n-6b) 

’ z = (v t ,e t + i ( e t .V t ) In tt 2 ) and 

(n-6c) 

h z = -i* V t .fit . 

(n-6d) 

p 


Substituting (II-6a) in (II-6b) we get 


— - J -( pVth z + kn 2 J~&zxV { J. (n- 7 ) 

JtV-/? 2 ' v / 


Substituting (ET-7) in (II-6a) we have 
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It can be seen from (II-7) and (II-8) that the- transverse fields are a function of 
the longitudinal fields only. Consequently if we can find the z-component of 
the fields, we automatically have the whole solution. 

II-3 The Homogeneous Vector Wave Equations 

In order to find the longitudinal components of the fields, we transform the 
Maxwell equations. Equations (II 3), into a more convenient form. We will 
carry out this derivation just for the H- field. The derivation for the E- field 

is similar. 

Using Equation (II 3b) we can write E in terms of H. Substituting this result 
in Equation (IT3a) we get 


V x 


V x H 


= k 2 H 


(II 9) 


For a step profile fiber with infinite cladding n is given by 


{ ttcore 0 < r < fl 

Nciad a < r < ~ 


m-io) 


where a is the core radius and r is the radial cylindrical coordinate. As it can 
be seen from (11-10), the index of refraction is constant throughout the fiber 
except at the core-cladding interface. Consequently, within the core and the 


II. Bound Modes 


23 


cladding. Equation (II 9) and its E- field counterpart will satisfy the equations 


V x (V x H ) = k n H 


(11-11 a) 


VxjVxE ) = knE 


OH lb) 


respectively. Expanding the left hand side of Equations (11-11) and substituting 
Equations (II-3c) and (II-3d) in the expansion, we get 


V H = -k 2 n 2 H 


(n-iic) 


V 2 E = -k z n 2 E 


(n-lld) 


We must emphasize that the operator in the left hand-side of the above 
- 2 

equations, V , is a vector operator, not to be confused with the scalar 

2 

Laplacian operator, V . In a curvilinear system of coordinates, the vector 
operator couples the components of the fields this makes 


- 2 2 
V *v . 


However, in the Cartesian system of coordinates. 


2 2 

V =v 
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[Snyder et. al., 1983]. We should warn that Jackson [1975] does not make this 

distinction in his Vector Formulas displayed inside the front cover. This 

2 2 

gives a false impression that the relationship V =V always holds. 

It must be stated that Equations (II 11) do not hold at the interface but only 
within the core and the cladding. If we want to know the field everywhere in 
the fiber, including the interface, we have to impose the boundary conditions 
of Maxwell’s equations in Equations (11-11) [Snyder et. al., 1983]. 

Equations (IM1) are known as the Homogeneous Vector Wave Equations 
(source-free). 

1 1-4 Solution of the Longitudinal Fields 


Since the fiber is cylindrical^ symmetric it is convenient to use a cylindrical 
coordinate system to solve Equations (11-11). The two vector differential 
equations obtained above can be broken down in six scalar differential 
equations. Four of them involve the transverse components of the fields 
which are coupled to each other in the equation. The other two involve the 
z-components of the fields only. As discussed before we can specify all fields 
in the fiber just by solving the longitudinal components. The equations that 
are of interest to us are expressed below 
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2 2 
d e z f l de z , 1 d e z {,2 2 n 2) 

~ T + 7Y ” - p l e *=° 

dr * r dip 


2 2 
d h z \ dh z 1 d h z | 2 2 2] 

+ — + \k n - P )h = 0 


■ + 


dr 


dr 


2 2 


r d* 


(II-12a) 


(II-12b) 


As it can be seen from above, the equation for h z has the same functional 


form of e z . For this reason we can solve our problem by treating only one of 


the equations. Defining new variables we can rewrite Equations (1112) as 


2 

d e. 
dR : 


2 

l de z 1 de z 2 

+ — 7 — r + u e z = °; 0<R<1 


R 


dR 


(III 3a) 


dip 


2 

d e 


dR ‘ 


2 l de z 1 3 e z 2 
+ - — + -r— --W e =0; kr< 


R 


dR 


R‘ 


(n-13b) 


dtp 


where R is the normalized radius R = r/a t 


U 


= a< 


k tt<?ore ft Slid 


(n-i3c) 


w = «V-tVi, d 


(II-13d) 


It is also useful to define 
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V = = fU 2 + W 2 . (H-13e) 

The V-number defined above is a measure of the number of modes that 
propagate in the fiber. The lower it is/ the lesser the number of propagating 
modes. On the other hand, the parameters U and W are a measure of how 
close to cut-off a particular mode is (we will see later on that each mode has a 
particular U and W values). Notice that their maximum value is V. A 
mode that is close to cut-off, has a U (W) value close to V (zero). In the 
same way, a mode that is further from cut-off, has a U (W) value close to 

zero (V). 

Equations (11-13) can be separated into the equations 

R 2 d 2 F , R dF + u 2 r 2 ,. l d 2< fi . = v 2 ; 0 <,R< 1 (II-14a) 

F dR 2 F dR d<p 2 


R 2 d 2 F , R d£ - VV 2 R 2 = - ^ 2< ^ - = v 2 ; 1 <R< 00 (II-14b) 

F dR 2 F dR <t> d(f> 2 

where we have made e z =F(R)<I>(0). By solving Equations (11-14) and 
imposing the condition of continuity of the longitudinal components across 
the interface, i.e., at R = 1 we finally find 


e z~ A v,n 


MVl 


fM 


h z =B v,n 


J v(^v,ji) 


8 v 


0 ^ R < 1 


(III 5a) 
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■'•■SI'-'* *>•«.»»»■*»« 


where 


cos 


(v^) 


/vW = 


l sin(v^) 


8 ^ 


-W= 


- si"M even modes 

cos (v$) °dd modes 


(H-15b) 


1 v (U v ,i) an d K V (W V ' are the Bessel function and the modified Bessel 
function respectively. 


The second subscripted index in A, U and W is consequence of further 
boundary conditions. We will see later that by imposing continuity of the 
tangential field components at the fiber interface we get an eigenvalue 
equation for each mode v. The fi index arises because of the possibility of 
multiple roots in the eigenvalue equation. The double index notation is also 
carried along to the propagation constant / J . 

II-5 Solution of the Transverse Fields 

Now we can determine the transverse components of the fields. They are 
found by substituting Equations (11-15) in (H-7) and (H-8). The final result is 
written below 
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Q Jv -1 (U v< ^ R)~ ^2 Jv + 1 (Uv,n R) f 
fv W 

U v , fi Q K v .,( VV v , R ) + C 2 Ky + 1 ( W v> p R ) 

v^7 Miy M ” 


0 ^ R < 1 


1 < R < 00 


£<jp 


Cl Jv - 1 (L/ V/ ,^C 2 /v + 1 (ii V/ /i ^ / \ 

/v (U v ,„) 8vW) 

Uv,n C] Ky _\{ W Vi n R)-C 2 K v+ ,{ Wy # n R) ^ i ^ 

w^7 * v 


0 < R < 1 

1 £ R < « 


hr = 






Qi /v - 1 (U v , n R) ■ Q Jv + 1 (Uv, /i R) / * 

Mu^j 

U v ,„ C 5 Kv.i(Wv^R) + C 6 Kv t l (Wv /M R) ^ /<|A 
Kv (W„) 


0 < R < 1 
1 ^ R < 00 


h* = 


CsJ v .\(U v ,n R) + C4j v + \(U V/ n R ) r 

1M fA(p) 

Uv , ± C 5 Kvi(W v , M R)-C 6 Ky + 1 (W v>p R) 
Wy <M Rv (W V/At ) 


0 < R £ 1 

1 < R < 00 


where 



(II- 16a) 


(H-16b) 


(II- 16c) 


(II-16d) 
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Q = 2aU~ ( ^ + ; 

Cs= 2M^[ !5v ^ Bv ^' Vf*^ d ^) ; 

Q = 2^7 Bv '^ + Vf*^ d • 

As discussed before, the components of the fields in the region R>1 is termed 
the evanescent field. 

The final solution reflects the existence of three kind of modes 

1) Transverse Electric Modes (TE Modes), v= 0 

2) Transverse Magnetic Modes (TM Modes), v = 0 

3) Hybrid Modes (EH and HE Modes), v* 0. 

TE Modes are modes whose e z components are null. TM Modes has a 
longitudinal magnetic field h z equal to zero. Hybrid Modes have both 
longitudinal electric and magnetic fields. They can be further classified in 
terms of relative amounts of E z and H z components [Snitzer, 1961]. 

II -6 Coefficients of Amplitude 

We can find the ratio of the coefficients of amplitude A Kfl and B v ^ by 
applying the boundary condition of the transverse fields 0at the interface. 
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Using the continuity of at R= 1 we get 


A v,fi 

B v ,n 


kui, w& net i fy-(u^) , i 

v/lv^V 2 ' S) W-J iIv(Ukii) Wvji K V (W V ,„)I 


(II-17a) 


Doing the same for the h p components we get 

B VJl kU v 2 „ W v 2 „ /arf n^,./;(»v^) + K >( tV -->-)| OT-17b) 

a k „ vP k „v 2 * n\Uv,„ !A u v) Wv.„ K v {w VJI )j 

The prime in the above equations represent a derivative with respect to the 
argument. As seen from above, we can either use 0117a) or (ll-17b) in our 
field solutions. As discussed by Snyder et. al. [19831, the amplitude A, „ could 
be set equal to the unity. However, it is chosen for consistency with the weak- 
guidance approximation (tt cor -tt clld ). Following Snyder’s normalization 
[Snyder et. al., 1983], we finally have 


i) for TE Modes 


€f — Cj 0 


_ 0£R < 1 

Ji(Uo,m) 

e *~\ K,(Wq, m R) 

\ K\(Wo.n) 


(II- 18a) 


1 < R < OO 


(n-18b) 
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0£R<, 1 
1 is R < °o 



0<>R<, 1 

I < R < oo 


h r =h z =e^=0 


e r = 


"L 


0 < R < 1 


1 < R < oo 


£ 


Z 


/oK„*) 

a Po,fi Jl( U 0.f) 

'"core ^o(W 0 /t R) 


0< K £ 1 


(II- 18c) 


(n-18d) 


(II-19a) 


(II 19b) 


1 ^ R < oo 


(II-19c) 
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0 < R < 1 


1 < R < oo 


iii) for Hybrid Modes 




U 

w 


KM 


KM 


”^ v .,(W v ^)-a 1 K^(W y)l R) i , , 


0<, R < 1 
1 <, R < oo 


a \ Jv-l(Uy,n R) ~ Jv +\(Uv,n g 0<R<1 

Jv {U v n) 6 

Uv,M <*\ K v-l(Wv,|i R) + fl2Kv + l(Wy, M R ) / x 1 ^ R < . 
W Vifl K v (Wv.„) 


iU v , n Jv (Ukm R) 


= 


a @KM J v ( Uv -^ 
iU KM (Wv,/i K) 

V fur pMW 

a Pv,n K V (W V(/1 ) 


/»W 0 < R ^ 1 


1 < R < 00 





£h fcH ArP fl 3 Tv - 1 (Uv, M R ) - fl 4 Jv ♦ 1 (jjVl jj R L (y) 

W) p v p Jv (Uv,,) 

Q kn&rs H3£ fl 5Kv.l(W y . M R)H-fl6Kv>l(Wv #M R) g (r) 

fJo p v ^ W V( , K v (IV v .„) 


(II-19d) 


(II-20a) 

(n-20b) 

(II-20c) 

0<R < 1 
1 < R < 00 


(II 20d) 
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/i«* = 


_ f knjh re /v - 1 {Uv, 

’ Pv,n J v {Uv'f*) 

£jg fc^re ^ v < £ g 5 - 1 ( Wy, p ^ ) ~ #6 Ky + 1 ( Wy, ^ R ) 




w, 




*v { W V ' n ) 


fvi<p) 
fv(<P) 


where 


K = 


i ^ Aa THs. a h I h iE^jj R ) „ 

v /« ka J v (U Kll ) 8v{ V] 


-I 


fig n F 2 K v { n R ) 


^v(<p) 


0 < R < 1 
1 < R < «x> 


, JF 2-1). „ _(F 2 + 1). 
ai JF 1-1). fll _(Fi + l). 

a , - F\_ + 1 ( Wclad f . 
5 2 2 '"core/ ' fl6 2 2 Wore) ' 

Uv.*W v.tf bi +("cU4j ncn Jb 2 ' f2 = 


V 


Uv.vWv.li) ^1 + ^2 
b _ 1 l Jy-i{Uy,n) Jv +i{Uy,n)\ 

2U VifJ \ Jy (U V 'n) Jy (U v ,fi) f 


and 


1 / K v - l ( W J , 

2V M K„( W vJ K v (w v j /. 


0 < R < 1 

1 < R < oo 

(II-20e) 

( 11-200 


Our problem will have solution only if the inverse of (II-17b) is equal to 
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(II-17a). This equality leads us to the eigenvalue equation and is discussed in 
the next section. 

II-7 The Eigenvalue Equations 

Right now we still do not know how to find the eigenvalues /3 VM , U Vfl and 
Yf v For the TM and TE modes they can be found in a straightforward way 
just by imposing continuity of the z-components of the fields at the interface 
(R=l). For the hybrid modes the eigenvalue equation is found from 
Equations (11-17). The resulting equations are transcendental in their 
arguments and are solved through numerical or graphical techniques. These 
equations, referred as the eigenvalue equations are displayed below for each 
mode [Snyder et. al., 1983; Snitzer, 1961]. 


i) TE Modes (continuity of h z ) 


/o(“oJ , w 
°^M ii 0 ,J °'" 


(II 21a) 


ii) TM Modes (continuity of e z ) 


k{UoJ +w ^ core =Q 


(II 21b) 


iii) Hybrid Modes (Equations (11-17)) 


II. Bound Modes 


35 


I Jy l U Kv) , ^cid K v'{ W Kh) \ 

\Uv,nJv {Uv,n) n&yre W Vrfl K v (W KM )J 

1 MHaA , K'(w k „) | / vA,„ \ 2 , v , 4 

^v,nJv (liv, /i) W v ,[iK v (W Kfl )l \kn COT e ) \U v ^W v ^j 


. (n-2ic) 


Equations (11-21) can be solved numerically with the help of Equations (II-13c), 
(II13d) and (II13e). 


It was pointed out by Snitzer [1961] that Equation (II-21c) is quadratic in 


KjUvud 

Uv'fijv (Uy r n) 


(n-22) 


Treating (11-22) as an independent variable. Equation (II 21c) can be broken 
down in two sets of solutions; one of them corresponds to the HE modes and 
the other to the EH modes. Their form is rather complicated and will not be 
displayed here. 


II-8 Weakly Guiding Fibers 

By definition, a weakly guiding fiber is a fiber which has an index of refraction 
of the core approximately equal to the index of refraction of the cladding 
[Snyder, 1983; Marcuse, 1974] 


Kore* *”dad “A 


(II 22) 
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Such fiber is almost similar to an unbound uniform medium. Being such it is 
a good approximation to assume that its fields are quasi Transverse 
Electromagnetic, quasi TEM, field. For this reason we use the quasi TEM 
approximation or 


e~h ~ 0 01-23) 

to solve the wave equation for the transverse fields. Using Equations (11-22) 
and (II-23) in the vector wave equation and solving it for the Cartesian 
components of the fields, we find that the transverse fields become solutions 
of the so called scalar wave equation. We want the results of the weakly 
guiding fiber only for posterior comparison with the exact solution. For this 
reason, we will limit ourselves to display them here [Snyder et. al., 1983] 


i) HE A+l M Modes (v=A+l, A2 0) 


| (x cos A <p - y sin XQ )Fa(R) ; 
£?tan \ (x sin X<f> + y cos A4 >)f*(R) ; 


even 

odd 


(II-24a) 


ii) Transverse Modes (v= 0, A=l) 


£tan “ 


(x cos <p + y sin <p )Fj (R) ; 
(x sin 0 - y cos 0 )Fi (R) ; 


TM 

TE 


(E-24b) 


iii) EH A+1/i Modes (v=A-l, A>1) 
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^tan 


I (x cos X<p + y sin X<p )Fa.(^) ; 
\ (xsinA0 -y cos X<p)Fx{R) ; 


even 

odd 


(n-24c) 


where 


FAR) = 


** K.„ r) . 

Ka(VVa,„) 


0<K£1 


l<R<oo 


The longitudinal z-components are almost zero and will not be displayed 
here. The eigenvalue equation for the weakly guiding fiber is 


,7. M M Ka.,(w A j ,) 

* h{^) •" X A (iv A ,) 


(11-25) 


where 


v = V uZ + ivj . 


The tilde in the eigenvalues was introduced in order to distinguish between 
the eigenvalues of the exact solution and the eigenvalues of the weakly 
guiding. The transverse magnetic fields can be found from 


h 


tan 



z xe 


tan 


(H-26) 


III. POWER INJECTION DUE 
TO SOURCES DISTRIBUTION 


III-l Introduction 

Fluorescent sources distributed in the cladding radiate incoherent light and 
can be modeled as radiating current dipoles with random phase and 
orientation [Christensen et. al., 1989; Snitzer, 1961; Marcuse, 1988; 1975]. As 
stated before, the total field injected into the core by these sources can be 
expressed as a sum over the bound modes. Each mode having an amplitude 

coefficient, c Vfl or 


E = 






( 01 - 1 ) 


Since the total field is excited by the current distribution j in the cladding, it 
is natural to expect that the amplitude coefficients, c Vfl are related to ;. 
Before deriving such a relationship, we will introduce a very important 
theorem: the reciprocity theorem. With the help of Maxwell’s Equations and 
the previous theorem, the amplitude coefficient follows [Snyder et. al., 1983]. 
Next we present a general equation for the power efficiency of a cylindrical 
fiber with sources distributed in its cladding [Marcuse, 1988; 1975]. Finally, we 
derive the power efficiency of a cylindrical optical fiber whose fields obey 
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Equations (1118) thru (11-20) and compare our results to the weakly guiding 
results of Marcuse [1988]. Our expressions are derived for the positively 
guiding fiber for two different distributions: 

1) sources uniformly distributed in the cladding and 

2) thin film sources distribution. 

For further applications of the reciprocity theorem, the reader should refer to 
Snyder et. al. [1983]. 

III-2 The Reciprocity Theorem (Conjugated Form) 

The conjugated form of the Reciprocity Theorem can be derived by defining a 
vector function F c 


F c = E xH * + E *x H 


(m-2) 


where * denotes complex conjugated. The unbarred fields satisfy Maxwell 
Equations with density current j or [Snyder et. al.,1983] 


V 


xE = i 



(III-3a) 


V xH=j n 2 E 
V fJo 


(ra-3b) 


The barred fields satisfy the conjugated form of those equations or 
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VxE‘ = H * 

V Go 

V xH * =/* + (« 'f 2 £ * 


(EII-4a) 


(in-4b) 


Using the identity for the divergence of two vectors we find 

V.F C = H *.(v x.E ) - E .(V x H *) -E*.(v xH‘)+H.(VxE‘) . (HI-5) 


Substituting (111-4) in (III 5) we get 


V.F c = -(e *. ; + E ((« »)2-n 2 )E.E» (HI-6) 


Applying the two-dimensional form of the divergence theorem to Equation 
(m-6) we have 


( V.F c dS =— f F c .zdS+ J> F c .ndl . (HI-7) 

Js dz Is Ji 

If we take S as the infinity cross-section of the cylindrical fibre, S^, the line 
integral then has to be taken along the circle r->«, where r is the cylindrical 
radius. Far from the core/cladding boundary, the amplitude of the fields of 
the bound modes fall-off exponentially consequently, at the line 

integral is null. Equation (III-7) reduces to [Snyder et. al., 1983] 
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I V.F C dS = — I F c z dS . (HI-8) 

Js. J s _ 

Equation (III-8) is the conjugated form of the reciprocity theorem. It relates the 
fields of two different fibers with indices of refraction n and n and current 
distribution j and j respectively. 

III-3 Excitation by Current Sources 

As it was stated before, the reciprocity theorem relates the fields of two 
different fibers. By knowing the fields of one fiber we may determine the 
fields in another. In our case, we would like to determine the forward 
propagating fields, E and H of a step profile fiber with sources distributed in 
it (Figure III-l). This fiber has a determined current distribution j and an 
index profile n. We choose a reference fiber with barred fields E and H, a 
current distribution j = 0 and an index profile W — n. Such reference fiber has 
a field solution given by Equations (IM8), (IM9) and (II 20) [Snyder et. al. 
1983]. Written in another way, we have 

E = e Ktl exp [i(3 v ,n z) ; H = h Vtfl exp (ip Vifl z) . (01-9) 

With the results above we can rewrite (01-6) and substitute it in (01 8). Doing 
that we obtain 
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Figure III-l. Sources in the cladding can be treated as current dipole vectors 
(arrows). Fluorescent molecules can be modeled as dipoles with random 
phase and orientation. 
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d_ 

dz 



F c .zd S = 


-1 


0 * 


■j exp (- i/J v/t z)dS . 


(m-io) 


The expansion coefficient in Equation (III- 1) is a function of z in the region 
z.,< z< z 2 and constant outside it. Substituting Equations (III-l) and (III-2) 
into Equation (III 10) and using the orthogonality condition 



x h a,p ) • zdS = 0 ; 


v * a 
or 

P *p 


(HI-11) 


we find 


dcv./i (z) 
dz 



exp z) dS 


(HI-12) 


where P Vfl is the normalization constant given by 


P 




x 


h v*n ) • z da . 


(HI-13) 


Finally integrating Equation (HI-12) we get 
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(III-14) 


which is a result obtained by Snyder et. al. [1983]. 

In any experiment, we should expect to collect the light in a region far away 
from the sources. Consequently, in the above result, we are specifically 
interested in the region z>z 2 . A slightly different expression can be obtained 
for backward propagating modes. 

III-4 Power Injected into Bound Modes 


The power carried by each mode can be expressed in terms of the amplitude 
and the normalization coefficients and written as 


jCv,/x| 2 Pvji 


(m-15) 
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Taking the sum over all modes of the previbus expression, we get the total 
power due to the bound modes, P core [Snyder et. al., 1983; Marcuse, 1988] 



V, n 


(HI-16) 


Since there are many sources distributed in the cladding and they have 
random phase and orientation, we can rewrite Equation (III 16) as 


(Pcore) — 



(DM 7) 


which is the ensemble average of the total power of the field [Marcuse, 1988], 

Equations (III-14) and (III-17) provide some clues on the characteristics of the 
amplitude of the injected modes. For instance, if j has only a z-component, 
the TE modes will not propagate. In the same way, if the z-component of the 
electric field of a particular mode is small, so will be the contribution of this 
mode to P core (the TM 0 mode of the anti-symmetric slab is a good example 
[Christensen et. al., 1989]). On the other hand, we could also choose the 
direction of j with respect to e Vfl in such a way to maximize P core However, 
a more realistic model for fluorescent sources is a random distribution where 
there is no preferred direction of orientation[Marcuse, 1988; Christensen et. 
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al., 1989]. In this case, we should expect every mode to be excited. 


As discussed before we should consider only that portion of the amplitude 
coefficient which is at the region z>z 2 . Rewriting it we have 


c v,n 




(m-i8) 


The ensemble average of the square of the modulus of the amplitude 
coefficient follows 

L ,4 = _L_/f e v .Jr)J( r)dv( <v> '). j • (r •) dV) . (m-19) 

Rearranging Equation (IH-19), we get 



.(/{r)y(r , )).«v J ,(r ')dVdv\ 


(HI-20) 


The quantity between parenthesis is a tensor of second rank. It can be 
represented by the matrix 



-4 -/ f [ < vV W -(j W ;• <r ')dV dv). (DI-21) 

lbPvji \l ' / 


Equation (in-20) can be further expressed as 
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e v,n ( r ) • (/ ( r ) j* (r )) (r ') dV dV 


(m-22) 


Equation (III-22) was previously obtained by Marcuse [1975]. The term 

(;• (r) i • (r •)) 


is the ensemble average of the current density tensor over many similar 
systems [Marcuse, 1988]. The current source is composed of many randomly 
phased and oriented dipoles. Being randomly phased and oriented, the 
components of j are uncorrelated. This makes the off-diagonal terms of the 
matrix null or 


/ 


[jx (r ) // (r ') ) 


0 


0 


(; (r) j * ( r ')) = 


0 {jy ( r ')) 0 

0 0 ih ( r ) jz* (r ')) 


\ 


/ 


(III-23a) 


In the same way, similar components of ; are also uncorrelated at different 
positions or, j { (r) and ;• (r r ) (i =x,y and z), are uncorrelated at r*r'. 
This makes the diagonal terms proportional to the Dirac delta function or 


(/» (r ) ji * ( r ')) « 1 ji ( r )| 2 S(r - r');i = x,y,z . 


(m-23b) 


Since there is no preferred direction of orientation, each component of j 
contribute the same amount to our result or 
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W=feHJ- < nI - 23c) 

This makes all diagonal terms in the matrix equal to each other. So we can 
assume Equation (IIT23a) to be proportional to the unitary tensor I . Thus we 
can write it as [Marcuse, 1975; 1988] 


(/' (r) (r ')) = S U(r-f'). 


(ffl-23d) 


Where S is the source strength. Substituting Equation (III 23d) into Equation 
(III-22) we get [Marcuse, 1975; 1988] 


b.,F) = -^rf s| e „„(rfdv 

lO'v./i Jv, ource 


(ni-24) 


which, in conjunction with Equation (III-17), can be used to find the power 
injected into bound modes 


P 


core “ 



v , \l 


(r f dV . 


(ni-25) 


Equation (III-24) can also be used to find the ensemble average of the 
amplitude coefficient of a continuum set of modes. 
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III-5 Total Power Radiated 

In this section we will derive an expression for the total power radiated by the 
sources in the cladding. The sources are assumed to be already excited. They 
are distributed uniformly in the cladding between the inner radius r in and 
the outer radius r out (Figure III-2) in such a way that a<r in <r out <b. The 
excited radiation is incoherent and its associated fields can be represented by 
an integral over a continuum of plane-wave modes. They can be expressed as 
[Marcuse, 1975; 1974] 



/ 2^P ( 

V JZ 


exp {-ik KtC .r)e K(T 
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where co is the angular frequency, P Ka is the mode normalization, e^is 
the unit vector that determine the direction of E K (Jf k K(T is the propagation 
vector and k, <t and /? are the x, y and z-components of the propagation 
vector respectively. 


The total power radiated is written as an integral over the continuum of 
modes [Marcuse, 1975; 1974]. Its form is similar to its discrete counterpart. 
Equation (IIT17), and can be expressed as 


Figure III-2. Sources uniformly distributed between radius r in and r 
(a<r in < r out <b). 
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where \c{K,of) is the ensemble average of the square of the modulus of the 
amplitude coefficient of the excited modes. It is given by [Marcuse, 1975] 


\c{K,of) 



E Ktr {rfdV. 


(HI-28) 


The factor of 2 in Equation (III-27) accounts for two possible polarizations of 
the transverse electromagnetic wave (TEM) [Jackson, 1975; Marcuse, 1975]. 
The z-component of the propagation constant is a function of the other 
components or [Marcuse, 1975] 


P = ^k 2 n 2 - k 2 - a 2 (HI-29) 

where k is the circular wave-number of the free-space excited radiation and 
n is the index of refraction of the cladding material. 

Substituting (III-26a) in (III-28) and assuming a cylindrical distribution of 
sources with inner radius r in and outer radius r we get 
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In the above formula, L, is the length of the fiber along which the fluorescent 
sources are distributed. 

We can find the total power radiated by substituting Equation (III-30) in 
(in-27) or 


Prad = 2— SL^(r| ut - 

32tt 2 




dfcdcr . 


(m-31) 


Equation (III-31) can be easily solved if we express the components of the 
propagation vector in a spherical coordinate system. In this system the 
components are expressed as 

K=nksindcos<p, o= n k sin0 sintp, p Ka =nkcosd (I3I-32a) 
dxrdf7=(« kf cost? sin# d#d0 . (III-32b) 

Substituting Equations (IH-32) into (m-31) we finally get [Marcuse, 1975; 1988] 

COfl o [ 2 2\ 

P rad = -^-”* SL \ r out- r inJ . (HI-33) 

Using Equations (III 33) and (III-25) and writing the angular frequency in 
terms of the propagation constant, we find the power efficiency 
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where the terms inside the sum are due to the bound modes. 


Notice that the above result is independent of the cladding radius b. This 
occurs because of the infinite cladding assumption used in the derivation of 
the fields in Chapter II. However, later on for the bulk distribution, we will 
make r out =b • Notice that this does not mean that the power efficiency is a 
function of the cladding radius b but a function of the outer radius of the 
source distribution. 


III-6 Power Efficiency 


For the exact case, the integral in Equation (III-34) is different for each mode. 
They involve the integration of the product of the normalized radius by the 
square of the Bessel K n function or RK n 2 (W v ^R) which is tabulated. The 
normalizations follow from Snyder et. al. [1983] whereas Equations (ni-35a), 
(III-36a) and (IH-37a) are the results of this work. 
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iii) EH and HE modes (v > 0) 
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(m-37b) 


where / is an integer and K out and R in are the normalized outer and inner 
radius of the fluorescent distribution given by r out /a and rj a respectively. 
They obey the inequality l£/? in <R 0Ut . Using the eigenvalue equation of the 
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exact solution of the TE and TM modes and Jhe recurrence relationship of the 
Bessel function, we can express the ratios ] V {U v ^) / J v+ 1 ( U v,/P 
J V (U Vfl )/J v _i(U v>fl ) in terms of K v < W v>fl ) / K v+ j( W Vfl ) and 
K v (W Vfl )/K v __ l (W VM ) respectively (see Marcuse [1974], page 68 for a similar 
derivation for the weakly guiding approximation). Substituting this result in 
Equations (III-35b) and (UI-36b) we get Equations (III-35c) and (III-36c) 
respectively. The equation of the power efficiency must take into account the 
degeneracy of the modes. Consequently, the final equation for the hybrid 
modes must be multiplied by two, odd and even modes. 

For the weakly guiding case the normalization and the integral of Equation 
(III-34) are the same for every mode. As before, the integral below can be 
determined in a way similar to the one above. Also, Equation (III-38a) is the 
result of this work. In the equations below we have used the expressions of 
the weakly guiding fields in Equations (11-24). 



— ^ — r (flout Ik^VVa,^ flout] ■ flin I^X. (W*,pR in]) 

AaIWa.^) 


(m-38a) 



where 
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I 1 for TM and TE modes 
A = \ v - 1 for HE modes 

1 v + 1 for EH modes 

We should stress the fact that the results of Equations (III-35a), (III-36a), (III- 
37a) and (III-38a) are valid for the general cylindrical source distribution 
depicted in Figure (III-2) assuming that the indices of refraction of the 
cladding and fluorescent sources are similar. Such assumption is necessary in 
order to avoid further complication in the expressions of the fields which 
arises whenever the fiber has more than two different indices of refraction 
(finite cladding). The fields solution for this fiber [Kuhn, 1974 and references 
therein] is helpful but will not be considered here. This assumption, that the 
fluorescent sources have to have an index of refraction similar to the cladding 
one, becomes unnecessary whenever the source distribution is either a thin 
film (r out =r in +5, S small) or a bulk distribution (r out =fc). In the first case, 
the index of refraction of the thin film can be neglected because of its small 
thickness. In the second case, the cladding is completely filled with the 
fluorescent sources which makes its index of refraction equal to the index of 
the fluorescent sources. 

Notice that the power efficiency of a weakly guiding fiber is dependent on 
only one index of refraction. As stated before. Equations (III-35) thru (III-38) 
are independent of the cladding radius, b, because of the infinite cladding 
assumption in the weakly guiding and exact solutions. Due to two distinct 
polarizations and the existence of odd and even modes. Equation (III-38a) 
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should be multiplied by two for A=0 and by four for A>0. 

Substituting A=1 (TE/TM modes) in Equation (III 38a) and having in mind 
that W 0 ^W li(1 when n cor( r"„ a d- we notice that Equation (III 35a) is 
already reduced to the weakly guiding case. Equations (III-35c), (III 36a) and 
(III-36c) can also be reduced by using the weakly guiding approximation 

P** ^ n cor e‘^ n c\ad 

For simplicity we will not write down the expression for the power efficiency. 

We will do so in the next sections. There, we will find the power efficiency for 
two particular distributions of sources: the bulk distribution and the thin 
film. Table (III-l) summarizes the results obtained in this Section. 

III-7 Power Efficiency of a Bulk Distribution of Sources 

An optical fiber with a bulk distribution of fluorescent sources is illustrated in 
Figure (I-3a). It consists of sources uniformly distributed in the cladding, from 
the core/cladding to the cladding/ third medium boundaries. It can be 
exemplified, but is not restricted to, a bare fiber core coated with a polymer 
which has fluorescent molecules dissolved in it, the so called optical fiber 
distributed sensor [Lieberman et. al., 1990; 1988; Blyler et. al., 1989]. If the 
radius of the cladding is big enough. Equations (m-35) thru (ffl-37) can be used 

to model it. 
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Choosing a fiber that has a cladding radius b, Equations (III-35) thru (HI-37) 
can be rewritten by making, r in =a and r out =&. This makes R in =1.0 and 
*out-*/«- ^ nice simplification arises when we make b very large, i.e, by 
substituting r oui =b-+ «*> and R out -»<» in Equations (III-35a), (III-36a), (III-37a) 
and (in-38a) we get 
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respectively. At infinite, the determinant 1 1^04)1 decays faster than R out . 
For this reason, all terms containing R out vanish in the limit R oul ->~. 


Substituting Equations (III-35b), (IH-36b), (ni-37b), (III-38b) and (IH-39) thru 
(111-42) into Equation (HI-34), we get the power efficiency 
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TABLE III-l 



MODES 

| M 2 dv 

JV source* 

ASSUMPTIONS 

Exact Solution 

TE 

TM 

Hybrid 

Equation (III-35a) 
Equation (III-36a) 
Equation (III-37a) 

A and B 

Weakly guiding 

TE, 

TM and 
Hybrid 

Equation (m-38a) 

A thru C 

Summary equations for the integral of the square of the electric field over the 
volume of a cylindrical distribution of sources. Inner and outer radial limits 
of the distribution are given by r in and r out (see Figure III-2). R in and R out 
are the normalized inner and outer radius and are given by r in /a and 

r Q Ja respectively. 





Assumptions: 


A. Infinite cladding radius in the fields solutions (&->«»). 

B. Index of refraction of the source material similar 
to index of refraction of the cladding (^ source”^ clad 

C. Weakly guiding approximation (« C ore ss ”clad^ 
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(III -43d). 


For simplicity, we have dropped the arguments of K, and J* • Although the 
power efficiency of the EH and HE modes have the same form, their 
eigenvalues are different. Equations (m-43a) thru (m-43c) are the result of this 
work. However, Equation (III 43d) is essentially similar to Formula (8) of 
Marcuse [1988]. As discussed before, we had to multiply the original result of 
Equation (III-43c) by two in order to account for both odd and even modes. 
Finally the total power efficiency is given by the sum 
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We should stress that Equations (in 43) base on at least two approximations: 


1) the infinite cladding approximation, which was introduced in the 

derivation of the fields in Chapter II and 

2) the infinite sources outer radius approximation introduced in 

Equations (III 39) thru (111-42). 

As discussed before, the first approximation results in expressions for the 
fields which are independent of the cladding radius b. It also results in a cut- 
off frequency equal to zero for the fundamental mode. Although the infinite 
cladding assumption is widely used, it is the finiteness of the cladding that 
imposes a non-zero cut-off frequency for the fundamental mode [Black et. al.. 
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1986 and references therein]. 

For a few modes, the second approximation leads to a peak-discontinuity near 
their cut-off frequency in the power efficiency diagram. Such behavior is 
more evident at low ^-numbers. However, at high V-numbers, the large 
number of modes tends to smooth out the whole distribution of power 
efficiency. As a consequence a sharp increase in the power efficiency near cut- 
off becomes less perceptible. Instead, by using Equations (III-35a), (III-36a) and 
(III 37a) we avoid this discontinuity in the power efficiency curve. Tables (III- 
2) and (III-3) summarizes the results of this Section. 

In conclusion, the power efficiency of a bulk distribution of sources can be 
expressed in the following functional form 

P eff = P eff^ r out =b ' "cor* "clad' *>• (01-45) 

We will show that the number of parameters required for the power 
efficiency of a thin film distribution will be reduced. 

III-8 Power Efficiency of a Thin Film Distribution of Sources 

An optical fiber with a thin film distribution of fluorescent sources is 
illustrated in Figure (I-3b). It consists of a thin layer of sources concentrated at 
a specific distance from the center of the fiber in the cladding region. It can be 
represented by, but is not restricted to, a bare fiber core coated with a thin film. 
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TABLE III-2 



MODES 

) M 2 dv 

* V too PCM 

ASSUMPTIONS 

Exact Solution 

TE 

TM 

Hybrid 

Equation (III-39) 
Equation (III-40) 
Equation (IH-41) 

A and B 

Weakly guiding 

TE, 

TM and 
Hybrid 

Equation (HI-42) 

A thru C 


Summary equations for the integral of the square of the electric field over the 
volume of a bulk distribution of sources. In this case, R in =10, r in =a and r out 
and R out -><». Also the source material is the cladding itself, i.e., 

^ source - ^ clad* 


Assumptions: 


A. Infinite cladding radius in the field solutions (&-»«>). 

B. Infinite outer radius of the source distribution 

G Weakly guiding approximation (n COT ~n clad ). 
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TABLE III-3 


MODES 

p cor/E rad 

ASSUMPTIONS 

Exact Solution 

TE 

TM 

Hybrid 

Equation (III-43a) 
Equation (III-43b) 
Equation (IH-43c) 

A and B 

Weakly guiding 

TE, 

TM and 
Hybrid 

Equation (ni-43d) 

A thru C 

Summary equations for the power efficiency of a bulk distribution of sources. 
The corresponding Equations are obtained by substituting Equations (in-39) 
thru (III-42) and Equations (III-35b), (III-36b), (III-37b) and (III-38b) into 
Equation (HI-34). In this case the source material is the cladding itself, i.e., 
^sourcf^ciad - Equation (IH-43d) is similar to Equation (8) of Marcuse [1988]. 


Assumptions: A. Infinite cladding radius in the field solutions (b—*<x>). 

B. Infinite outer radius of the source distribution 
^out -400 ) in Equations (III-39) thru (HI-42). 

C. Weakly guiding approximation (» core -« clad ). 
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In this case, the air itself acts as an infinite radius cladding. 

The power efficiency of a thin film can be obtained by writing r out in terms of 
r where r L -r = 8 is the layer thickness of the film. For a thin film, the 
thickness 8 has to be much lesser than a/VJ v ^. Substituting 
R out = R in +8/a into Equations (III-35a), (III 36a), (III-37a) and (III-38a) we can 
find their corresponding thin-film expressions. These equations are a 
function of the determinant of K,(AR 0Ut ) which can be found by expanding 
it in terms of powers of 8. The first order approximation for the above 
determinant is given by 

IK, (A Rou.l-IK. (A Hint + AS/ a d (n*- 46 ) 

Again, using the first order approximation in 8 for R out and substituting 
Equation (III-46) into Equations (III-35a), (m-36a), (in-37a) and (III-38a) for the 
thin film case, we get the results below 
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which are also a first order approximation in 8. 


Neglecting the second higher order terms in 8 of Equation (ffl-33) we get 
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Notice that Equations (III-47) and (III-48) are all linear in 5. Consequently, the 
result of P eff is independent of this parameter. If we substitute these 
expressions into Equation (III-34), we get the power efficiency of a thin film 
located anywhere in the cladding. However, of particular interest, is the 
optical fiber which has a thin film located in the core/cladding boundary, i.e., 
£. n= 1.0. It is easier to make such fiber and it is also the one that has the 
highest power efficiency. Equations (III-47) and (III-48) can be easily reduced to 

this case. 

Using Equations (III-48) and (III-25) and substituting R in =10 in it, we get the 
power efficiency of a thin film distribution at the core/cladding interface 



Although not published, Equation (III-49) was obtained previously by 
Marcuse. 


The expresssion 
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can be simplified if we expand it and use the recurrence relationship to 
eliminate the highest and lowest order terms of the Bessel function. The 
arithmethic is very involved. The final result is 


d ~^ Tw~ + 2 1* = 2 k ' 2 C • 


(m-50) 


Substituting Equation (III-50) into Equations (III-47) we get 
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iii) EH and HE modes (v > 0) 
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iv) weakly guiding 
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Remember that, for the thin-film, the S in the equations above is very small. 
Finally, substituting Equations (HI-51), (III 35c), (III 36b), (nF37b) and (in 38b) 
into Equation (ni-49) and taking the limit, 5->0, we get the power efficiency 
of a thin film distribution over the core /cladding boundary or 


core 


rad 


i 


ITH 2 n dad fc fl V 


X 


u o., K . ? ( w °.J 


(ni-52a) 


core 


rad 


i n 


core 


i™ 


X 


P 0,a 


1 + 


a 2 t>l j 


Ji(Uq.m) 1 "L| k .( W <J 
/, 2 (U J »L *>„.„) I 


(ffl- 52b) 


core 


rad 


EH 2 ,3 3 

HE n clad n co^ a 


X 


III. Power Injection 71 


I 





2 + 2 + , 2 



W W a 2 ft 





K „( W J 


fl l a 3 


Jy -1 + ^^4) Jy -h l[ ^y,/j ^l^S^y - 1| ^y + l| 


U U J 


w. 




M W J 


<m-52c) 


core 


rad 


1 


weakly 2 2 2 2 

guiding ® ^ clad * * 


I 


V'H 


(2 - J 



k a-.I 

'w 


ik^.i 

K4 


For simplicity, we have dropped the arguments of K t and J t . The total power 
efficiency is given by Equation (in-44). 

Equations (ffl-52a) thru (ni-52c) are the result of this work whereas the weakly 
guiding case. Equation (ffl-52d), is similar to that derived by Marcuse [1988]. 

As it was summarized in Tables (III-4) thru (III-6) all Equations obtained 
involve at least one approximation -the infinite cladding radius 
approximation of the field solution. However, for the special case in which 
the cladding of our thin-film fiber is the outside air, « C ] ac j=l-0, Equations (III- 
52a) thru (IH-52c) involve no approximation at all. That is because the outside 
medium, air, essentially extends to infinity. The power efficiency P eff of a 
thin film coated fiber is a function of 
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P eff = r in' "corJ ”clad' *>• (III ’ 53) 

For a bare core fiber coated in the core/air interface r in =a and n clad =1.0 and 
the number of variables is reduced by two. Notice that in both the case of a 
thin film and a bulk distribution of sources the power efficiency is a function 
of many variables. This complex functionality makes it difficult to analyze the 
behavior of P eff . It is also difficult to sketch general predictions for a high 
injection efficiency by a fluorescent cladding optical fiber. In the next chapter, 
some results of model calculations based upon the above expressions allowed 
us to formulate a few rules for designing an optical fiber sensor. 
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TABLE III-4 


MODES 



M 2 dV 


ASSUMPTIONS 


TE Equation (III-47a) 

Exact Solution TM Equation (III-47b) A 

Hybrid Equation (HI-47c) 


TE, 

Weakly guiding TM and Equation (ffl-47d) A and B 

Hybrid 


Summary equations for the integral of the square of the electric field of a thin 
film source at an arbitrary position in the cladding. Inner and outer radial 
limits of the distribution are given by r in and r out = r in +<5 (<5->0). R in and 
K out are norma hzed inner and outer radius and are given by r in /a and 
r 0 Ja respectively. 


Assumptions: A. Infinite cladding radius in the field solutions (b—>°°) 4 

B. Weakly guiding approximation (n core =n clad ). 
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TABLE III-5 



MODES 

f I*?,,/ dV ASSUMPTIONS 

Exact Solution 

TE 

TM 

Hybrid 

Equation (III-51a) 
Equation (III-51b) A 

Equation (HI-51c) 

Weakly guiding 

TE, TM 

and 

Hybrid 

Equation (III-51d) A and B 

Summary equations for the integral of the square of the electric field of a thin 
film source at the core /cladding boundary. Inner and outer radial limits of the 
distribution are given by r in =<J and r out = a+5 (5->0). R in and R out are the 
normalized inner and outer radius and are given by 1.0 and 1.0+8/ a 
respectively. 

Assumptions: 

A. Infinite cladding radius in the field solutions (b— »«*>). 

B. Weakly guiding approximation ( « core =n c iad )• 
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TABLE III-6 



MODES 

P co,/P„d 

ASSUMPTIONS 

Exact Solution 

TE 

TM 

Hybrid 

Equation (III-52a) 
Equation (HI-52b) 
Equation (DI-52c) 

A 

Weakly guiding 

TE, 

TM and 
Hybrid 

Equation (IH-52d) 

A and B 

Summary equations for the power efficiency of a thin film distribution of 


sources. The corresponding Equations are obtained by substituting Equations 
(HI-51), (m-35c), (III-36b), (IH-37b) and (HI-38b) into Equation (in-49). Equation 
(HI-52d) is similar to Equation (8) of Marcuse [1988], 


Assumptions: 


A. Infinite cladding radius in the field solutions (M»). 

B. Weakly guiding approximation (w core ~« c i ad ). 


IV. POWER EFFICIENCY OF A 
DISTRIBUTION OF SOURCES 


IV-1 Introduction 

In this chapter we have used the results of the power efficiency of the 
previous chapter to determine its behavior as a function of different 
parameters and their combinations for specific conditions. 

As it can be seem from Equations (IH- 43a) thru (III 43c) and Equations (III-52a) 
thru (III-52c) the functional behavior of the power efficiency is too involved 
to be determined by direct analysis. Added to this, there is also the need to 
determine the eigenvalues of our problem from a complicated transcendental 
equation. A problem of this magnitude requires the use of computational 
techniques. For this reason, we have written a program in FORTRAN to 
compute the P eff against many of the possible variables. In certain cases, the 
parametric values used to determine the power efficiency are not realistic. 
This approach was used in order to determine the general behavior of the 
power efficiency. However, some of the parametric values are within the 

realm of the real world. 

We have used Equations (III-43a) thru (III-43c) for the bulk and Equations 
(IIT35) thru (III-37) for both bulk and thin film distributions. As discussed 
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before, if we use Equations (III 35) thru (III ^7) into Equations (III-34)we get a 
more exact and general result. I.e., by using an appropriate R out and R jn we 
can determine the P eff for any distribution of sources with cylindrical 
symmetry. For the bulk distribution, the difference between the two 
approaches lies in the upper radial limit used in the integral of the square of 
the electric field of the bulk distribution. The first one used an infinite value 
for the upper radial limit or 


Jt 


vJt dv = 


III 


\ e v,$ rdrd<pdz 


=r k> 


and the second one used the actual value of the outer radius of this 
distribution or 


1 


v,#dV 


■in: 


\e v ,$ rdrdcpdz 


These calculations were performed differently in a heuristic approach. 


The eigenvalues U V fl and W Vtfl are found by using four different 
transcendental equations for TM, TE, HE and EH modes. They have limits 
that are either well-known or can be calculated [Snyder et. al., 1983]. Using 
these limits, we can square the eigenvalues. The ZBRENT subroutine of 
Numerical Recipes [Press et. al., 1986] was used to find them. This subroutine 
employes the Secant method in conjunction with the Bisection method. 
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As it was noticed before, the power efficiency of both bulk and thin film 
distributions are very complicated. It is a function of up to five parameters 
which involve different combinations of fl, r out , ti cortf ti c j ac j and X. In 
order to predict its behavior, we have used two different approaches: in the 
first one we have analyzed the behavior of the P eff for a variable V-number 
and in the second one for a constant V-number. 


In the first approach, the V-number 

V = a fcV Me?, re - N c ?ad- tli 2 + W 2 

can be varied by changing one of its parameters separately namely, a, X, 

”core and n ciad- The results were P lotted either a g ainst the ^-number or 
against the actual variable. All of them have revealed a jagged appearance 
more or less perceptible depending on whether we are treating the thin-film 
or the bulk distribution (Figures (IV-1) thru (IV-4) and Figures (IV-6) thru 
(IV-11)). Marcuse [1988] was the first one to discuss this sharp increase. In his 
paper he has plotted the power efficiency of a weakly guiding fiber against the 
V-number for a fixed core/cladding radius and wavelength. He has also 
pointed out that this sharp increase occurs near the cut-off values of the 
modes. Near cut-off, the evanescent field is stronger and there is a higher 
probability of power injection. The weakly guiding formulas obtained is 
simple enough to allow him to vary the V-number freely, with no regard to 
the other parameters. Although, in the real world, his graphs are equivalent 
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to a plot of the power efficiency against the index/indices of refraction, he has 
not treated neither this case nor the remaining ones, i.e., variable wavelength 

and core radius. All these cases are treated in this work and the results are 
displayed in this chapter. 

In the second approach the behavior of the P eff was analyzed for a constant 
V-number. The V-number was kept constant by changing only two of its 
parameters. It was used two techniques: 

1) to vary the core and cladding radius a and b and the wavelength A 
at a constant ratios a / A and b/ a and constant indices of refraction n 

core 

and ”dad (notice that although we are varying three parameters only two are 
a V-number parameter). It is interesting to notice that for this specific case, 
the power efficiency remained constant for variable wavelength and core 
radius. This result suggests that both a / A and b/a are two independent 
variables. Notice that the expressions of the thin-film do not involve the 
cladding radius b. Consequently, a similar conclusion which does not 
involve b and combinations of it also applies. This result essentially 
decreases the dimension of the problem by one variable. This makes its 
analysis simpler and 

2) to vary the indices of refraction of the core and cladding, n core and 
n ciad respectively, at the same time. Remember that there are specific values 
for the indices of refraction « core and n c]ad that allow a constant V-value. 
The final result was plotted against the difference M core -n clad and revealed an 
almost linear increase in a log-log scale. 
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The constant V approach can be justified by the fact that for a constant V- 
number the power efficiency curve is smoother (there are no sharp peaks). 
This makes it easier to fit a polynomial equation to the graphs and analyze its 
behavior. Notice that since the eigenvalues of a weakly guiding fiber are 
independent of the fiber parameters, the P eff of this fiber at constant V- 
number is directly proportional to the coefficient of the terms inside the 
summation sign (Equations (III-43d) and (III-52d)). However, this 
simplification does not occur in the expressions of P eff derived in this work 
(Equations (m-43a) thru (III-43c) and Equations (III-52a) thru (III-52c)). This is 
because the eigenvalue equation of the exact solution is more complicated 
and dependent on more parameters. 


IV-2 Power Efficiency at Constant a/X 

By varying a and X in such a way to make the ratio a/X constant and 
fixing the indices of refraction, n core and n clad and the normalized- radius, 
R t and R in (b/a and 1.0 respectively for the bulk distribution), it can be 
concluded that the P eff remains unchanged. In order to see that, first we 
should recall the values of U V fi , W vfl and V or 


U 




= a 


k ^H(?ore Pv,n ' 


(n-i3c) 


w 


/ t n — a Pv,n ~ k tt c ^ ad and 


(II-13d) 
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V = akj K c ?> r e - H c ? ad = 1 U 2 + W 2 . (II-13e) 

Notice that /? can be rewritten as 

Pv,n = fc^efc (IV-1) 

where n efc is the effective index of refraction which obeys the inequality 


^clad — ^efc ^ ^core • 


Since A: = 2k/ X, it can be easily seen that U V fl , W v and V are the same 
for the prescribed conditions 


a - a \ 
X A, 


- £n. 
An 


(IV-2a) 


and 


”core' M ciad ' K in and R out constant. (IV-2b) 

Finally, substituting Equation (IV-1) into Equations (ffl-43) and (III 52) under 
Conditions (IV-2), writing b as a proportion in a and having in mind the 
previous result for U Vfl , W Vfi and V, we find that the P eff is the same for 
different values of a and X. The same conclusion also applies whenever we 
use Equations (HI-35a), (ni-36a) and (IIT37a) in Equations (HI-34) of the power 
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efficiency using the relations or r out = fl ^out an< ^ r in~ consequence 

the functional dependence of the power efficiency of both bulk and thin-film 
distribution are reduced to 

P e ff = P ^core' ^clad^ - 
P e ff = ^in' n core' ^clad^ 

respectively. For the special case of a thin film distributed in the core/cladding 
boundary, Pj n =l 0, the previous functionality dependence is reduced to 

P e ff = n core' n clad^ 

IV-3 Power Efficiency of a Weakly Guiding Fiber 

The results for the weakly guiding case are presented in Figures (IV-1) thru 
(IV-4). They were derived using the same values as those in Marcuse [1988], 
namely, >1=1.3 pm, «= lO-O^m and b/a= 5.0. In order to vary V, we used 
a fixed index of refraction for the core and varied the index of refraction of the 
cladding, Hdad. The difference ep=n core -n clad had to be small enough to 
simulate a weakly guiding fiber (ep< 0.1). For Figures (IV-la), (IV-3) and 
(IV-4a) ep varied roughly from 10' 3 at V= 0.05 to 10‘ 2 at V= 10. For Figure 
(IV-2) it varied from 10' 7 to 10' 2 . For higher V-numbers ep would increase 
even more invalidating the weak guidance condition. For this reason, we 
have avoided going further into a higher V-number region. For sake of 
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comparision, we have reproduced with .permission Figures 2 and 3 of 
Marcuse [1988] (Figures (IV-lb) and (IV-4b), respectively). The figures are for 
the bulk distribution. Figures (IV 1) thru (IV-3), and thin film. Figures (IV-4). 

i) Bulk Distribution 

We have plotted three graphs for the bulk distribution. Figures (IV-la) and 
(IV-2) follow from Equations (III 43a) thru (III 43c). Figure (IV-3) was obtained 
using the actual radial outer limit b of the source distribution (Equations 
III 35 thru m-37). We start by discussing Figures (IV 1) and (IV-2). 

The labels of Figures (IV-la) and (IV-2) correspond to the V- values at which 
the modes start to propagate. For instance, in a weakly-guiding fiber the 
modesTE 01 , TM Q x and HE 2 ^ have a cut-off at V =2.405. This is the 
corresponding V-value for the first sharp peak in Figures (IV-1) and (IV-2). 
All plotted peaks occur at cut-off of the labeled modes [Marcuse, 1988]. 
Between 0 and 2.405 only the fundamental mode, HE 1 j, propagates, 
consequently only this mode contributes to the power efficiency in this 
region. The plateau observed between V=0 and 0.5 occurs because between 
these values 


U^o,\ = V. (IV-3) 

The height of each plateau can be found by substituting Equation (IV-3) into 
(III-43d) or 
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Figure IV-la. Power Efficiency of a bulk distribution of fluorescent sources 
versus the V-number (weakly guiding case). The wavelength, the index of 
refraction of the core and the core and cladding radius are held fixed at 
A=1.3jum, n COT =\A, a = lOjim and b=50.0^m respectively. The Index of 
refraction of the cladding varied from 1.399 at P=0.05 to 1.385 at 1/=10.0. 
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Figure IV-lb. Power Efficiency of a bulk distribution of fluorescent sources 
versus the V-number using the weakly guiding approximation. The 
wavelength and the core and cladding radii are held fixed at A=1.3^m, 
a=\0.0jim and b= 50.0/im respectively. (From D. Marcuse, "Launching 
Light into Fiber Cores from Sources Located in the Cladding", IEEE 
Journal of Lightwave Technology , Vol. 6, No. 8, Aug. 1988, pp. 1273, 
reproduced with permission; © 1988 IEEE). 
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P core 
P rad 


0 <V <05 
A^0,/i = 1 


1 

(b 2 - a 1 ) « c 2 ,ad 


\ Kq 2 (W w ) 

K? ( Whi ) 


1 . (IV -4) 

[b 1 - a 2 ) n c 2 lad fc 2 


where 


Wa = o,i = 0, 


lim Ko (z) = - In z + In 2 - 0.577 

z — > 0 


and 


z — > o 2 

For « clad =l-4 Equation (IV-4) becomes 9-lxlO -6 . 

Notice that there is a cut-off degeneracy among many modes in a weakly 
guiding fiber. For instance, the modes HE 2 l fl , TEq M and TMq >fi have the same 
cut-off frequency. However, for an arbitrary fiber, part of this degeneracy is 

removed. 

The results in Figure (IV-la) are essentially the same as those of. Figure (IV- 
lb) [Marcuse, 1988]. However, the results of this work are slightly higher. We 

attribute this difference, possibly, to three factors: 
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1) we may have used a higher density of points than Marcuse did. By 
doing so we calculated more points near cut-off. As discussed by the previous 
author, this would give us a higher value for the power efficiency near the 
cut-off region (the closer to cut-off the higher the power efficiency), 

2) Marcuse's results assumed the weakly guiding approximation 
”core~ n ciad~” over the whole interval of V-values. For this reason, his 
result is independent of the two indices of refraction and, in order to obtain 
his figure, all he had to do was to vary the V-value only. In our case, in 
order to vary V we had to vary either n COTe or w clad or both. This results in 
bigger differences ep for bigger V-values. This difference, makes the final 
result of the power efficiency bigger than using the weakly guiding 
approximation and 

3) the index of refraction used in Marcuse [1988] might have been slightly 
higher than the one used in this work. This conclusion can be justified by 
analysing Equation (III-43d). Within the weak guidance approximation U 

V t /Z 

and W v> ^ depend only on V. So, for a fixed V-number, the sum in 
Equation (III-43-d) is always equal no matter what are the values of the other 
parameters of the weakly guiding fiber [Marcuse, 1988]. Consequently, at a 
specific V-number, the power efficiency of a weakly guiding fiber would be 
directly proportional to the square of the wavelength and inversely 
proportional to the square of the index of refraction and the difference of the 
square of the cladding and core radius. Stated in another way 
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P rad 


weakly 

guiding 


(i> : 


i 2 )"/ 


lad 


(IV-5) 


Equation (IV-5) does not apply when b=fl. 

Figure (IV-2) was obtained for n core =2 .6. The plot was made between V= 0 
to V= 10 so we could compare it more easily with Figure (IV-la). It can easily 
be seen that the power efficiency of the first fiber is 3.45 times greater than the 
power efficiency of the second. This result is in agreement with Equation 
(IV-5) which gives us the ratio 


=( 2 - 6 /l.4f = 345 . 


Figure (IV-3) was plotted so we could analyze the effect of the infinite outer 
radius approximation introduced in the integral of the square of the electric 
field. We can see from there that the peaks corresponding to the HE t n modes 
are slightly lower. Also notice that the plateau observed in the previous 
figures does not appear here suggesting that this feature is an artifact 
introduced by the previous approximation. A closer analysis of the numerical 
data indicates that, for some reason, this difference is restricted only to the 
cut-off region of the HEj n modes. In conclusion, the infinite outer radius 
approximation influences only the region near cut-off for the HE lrt modes 
and does not introduce serious errors in the results of this work. 
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ii) Thin Film Distribution 

The result for the thin film distribution is presented in Figure (IV-4a). In this 
graph, we have used a index of refraction « clad =l-46. The thin film is 
assumed to be located in the core/cladding interface. Compared with the bulk 
case, it is much smoother. The reason for this smoothness was discussed by 
Marcuse [1988] and is related to the evanescent wave that is emitted from the 
sources located further from the core/cladding interface. In other words, in a 
bulk distribution we have several layers of fluorescent material which are at 
increasing distances from the core boundary. The layers contribute more 
power near cut-off which produces the sharp increase in power in these 
regions. Therefore, in a bulk distribution, the power efficiency near cut-off is 
the sum of the contributions of many infinitesimal layers. This yields the 
discontinuity observed. However, for the thin film, we have only one layer. 
Consequently the integrated contribution near cut-off is not as dramatic. 
However we still have small peaks that occur nearby cut-off. 

It is clear from the figures that a thin film distribution in the core/cladding 
boundary is more efficient than a bulk one. However the total power 
radiated by the thin film is lower because there are fewer sources. In addition 
to that, the nearer the sources are to the core/cladding interface the more 
evanescent injection we get because the evanescent wave intensity decays 
exponentially from the sources. In this way, the closer the sources to the 
core/ cladding interface, the higher the probability of exciting bound modes. 
This also results in a higher efficiency. Our graph reproduces very well 
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Figure IV-2. Power Efficiency of a bulk distribution of fluorescent sources 
versus the V-number (weakly guiding case). The wavelength, the index of 
refraction of the core and the core and cladding radius are held fixed at 
A-1.3/im, n core = 2.6, a- 10/im and b=50.0/im respectively. The Index of 
refraction of the cladding varied from 2.5999 at V = 0.05 to 2.592 at V = 
10 . 0 . 
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Marcuse's results. Figure (IV-4b), at low V-number [Marcuse, 1988]. 
However, at high l/-number, it starts to diverge toward higher values for 
reasons discussed previously. 

In summary, we have used the exact field solution of the step index profile 
fiber to calculate the power injection efficiency of a bulk distribution (Figures 
(IV-la), (IV-2) and (IV-3)) and thin film distribution (Figure (IV-4a)). Figures 
(IV-la) and (IV-2) were obtained using Equations (III-43a) thru (III-43c). 
However, Figure (IV-3) was obtained substituting Equations (III-35a), (III-36a) 
and (III-37a) in the general equation of the power efficiency (Equation (HI-34)). 
We stress the fact that although Figures (IV-la) and (IV-3) are very similar, 
they were obtained using two different upper radial limits, r out — and 
r out = ^ respectively, in the integral of the square of the electric field. The 
reason for the previous approximation resides in the simplification 
introduced in the final equations without introducing serious errors. Figure 
(IV-4a) was obtained using Equations (III-52a) thru (III-52c). All Figures agree 
very well with the weakly guiding results obtained previously by Marcuse 
[1988] (see Figures (IV-lb) and (IV-4b)). He has obtained his bulk distribution 
results (Equation (III-43d)), using the infinite outer radius limit, r out — ■ ► in 

the volume integral of the square of the electric field (Equation III-42). Table 
(IV-1) summarizes the approximations involved in this Section. 
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Figure IV-3. Power Efficiency of a bulk distribution of fluorescent sources 
versus the V-number (weakly guiding case). The wavelength, the index of 
refraction of the core and the core and cladding radius are held fixed at 
A=1.3/im, n core =1.4, a = 10/im and b=50.0jxm respectively. The Index of 
refraction of the cladding varied from 1.399 at V= 0.05 to 1.385 at V= 10.0. 
This graph was obtained using a finite value, b, for the outer radius of the 
distribution in the equation for the power efficiency. 
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Figure IV-4a. Power efficiency of a thin film distribution of sources in the 
core/cladding boundary versus the V-number (weakly guiding case). The 
wavelength, the index of refraction of the core and the core radius are held 
fixed at A=1.3/xm, « core =l-46, and a=10//m respectively. The index of 
refraction of the cladding varied from 1.4599 at V= 0.05 to 1.322 at V= 
29.95. 
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Figure IV-4b. Power Efficiency of a thin film distribution of fluorescent 
sources versus the V-number using the weakly guiding approximation. 
The wavelength and the core and cladding radii are held fixed at 
A=1.3^m, a=10.0^m and b=50.0/xm respectively. (From D. Marcuse, 
"Launching Light into Fiber Cores from Sources Located in the Cladding", 
IEEE Journal of Lightwave Technology , Vol. 6, No. 8, Aug. 1988, pp. 1273, 
reproduced with permission; © 1988 IEEE). 
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IV-4 Thin Film at an Arbitrary Distance 


In order to determine the best position of the thin film within the cladding, 
we have plotted the logarithm of P eff against the normalized position of the 
film R in = r in /a (Figure (IV 5)) for a= 5.0pm, A=0.4pm, n cor =1.6 and 
”ciad =1 - 4 - We can see that the P eff decays exponentially with the distance 
from the boundary. Such behavior is expected since the total power associated 
with the evanescent field also decays in an exponential fashion. In 
consequence the closer the cladding sources are to the core/cladding boundary 
the higher the power efficiency. 

IV-5 Power Efficiency Versus Wavelength 

We have computed the power efficiency for TE 0 TM 0 ^ EHy and 
^i,/i m °des (p =1, 2, ... and v =1, 2, 3, ...) for a wide range of wavelength. 
The results are presented in Figures (IV-6) thru (IV-9). We have plotted the 
same data first as a function of the V-number and then as a function of the 
wavelength A. Notice that the actual independent variable is the 
wavelength. All other parameters are held constant. Since the V-number is 
inversely proportional to A, the left-hand side of each graph which have V- 
number as independent variable, corresponds to the highest wavelength, 
A=2.0pm. Similarly, the right-hand side corresponds to the lowest 
wavelength and vice-versa. As before, we present results for a bulk 
distribution and thin film. 
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TABLE IV-1 


Figure (IV-lb) Figure (IV-1 a) Figure (IV-3) 


A A A 

B C C 

D D E 


Comparative table on the assumptions used in three different results for the 
bulk distribution. 


Assumptions: A. Infinite cladding radius in the fields solutions (b-*~). 

B. Weakly guiding approximation (« core =n clad ). 

C. Arbitrary indices of refraction (n COT( ?n c]ad ). 

D. Infinite outer radius of the source distribution in the 
integral of the square of the electric field (r out ->°°)- 

E. Finite outer radius of the source distribution in the 
integral of the square of the electric field (r out = fr)- 
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i) Bulk Distribution 

We have plotted four graphs for the bulk distribution. Figures (IV-6) and (IV- 
7). In all of them, we have used r jn =a=5.0//m and r out =b=25.0//m for the 
inner and outer radius of the distribution respectively. However, we have 
used tt core =1.6 and n clad =1.55 in Figures (IV-6) we used and n core =3.0 and 
n dad = l-l Figures (IV-7). We have found that 

1) for a given wavelength X, the bigger the difference between « core and 
n c)ad the higher the power efficiency (we will show later that the same holds 
true for the V -number). It seems that a similar result (that the gain 
coefficient increases with the difference between the indices of refraction) was 
obtained for the TE 0//1 modes by Watanabe et. al. (see Figure 5 [Watanabe et. 
al., 1973]). This result can be easily explained for the TMq ^ modes in terms of 
the amplitude of the electric field. The electric field of TM 0 ^ modes in the 
cladding region is directly proportional to the square of the ratio n COT Jn c , ad 
[Snyder et. al., 1983]. Since the expansion coefficients are directly proportional 
to the amplitude of the electric field so will be the power injected into the core 
(see Equation (III-17)). No attempt was made to see whether there is a 
similarrelationship between the amplitude of the electric field of the hybrid 
modes and the ratio n cor{ /n dad - 
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Figure IV-5. Logarithm of the power efficiency of a thin-film distribution of 
sources versus the thin film normalized distance, R jn = r in / a. The 
wavelength, the core radius and the index of refraction of the core are held 
fixed at A=0.4/tm, a = 5. 0/im, n core =l-6 and « c i ad =l-4. The power 
efficiency is higher when the thin film is at the core/ cladding boundary. 
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2) There is a reinforcement in the power efficiency among the TE 0 
™o,m anc * modes of the weakly guiding fiber (« core =rt clad ) and a 

reinforcement between and HE 2 M modes of the strongly guiding fiber 

^cor^ ,>>tt ciad^- Reinforcement occurs because the corresponding modes 
have the same cut-off frequency in each guidance limit (see the cut-off 
equation of the hybrid modes in Table 12-4(b) of Snyder et. al.). In other words 
the corresponding modes have a degenerate cut-off frequency which can be 
removed whenever we use arbitrary indices of refraction, « cor( p« clad . Since 

they have the same cut-off frequency the power efficiency is necessarily 
reinforced, 

3) In general the longer the wavelength the higher the power efficiency. 
Apparently, this result reflects the characteristics of the behavior of a wave. In 
other words the longer the wavelength with respect to the dimensions of the 
fiber the more tunneling one should expect from the cladding. This also 
implies that the lower the V- number the higher the power efficiency, a 
result which is contrary to the previous belief that a higher V-number 
would yield a higher power efficiency [Marcuse, 1988; Lieberman et. al., 1990]. 
Although, the curve of the power efficiency in Figures (IV-1) thru (IV-3) does 
increase with the V-number, it must be stated that in those cases, V was 
obtained by changing the difference between the indices of refraction while 
here we have changed the V-number by changing the wavelength. In 
conclusion, the behavior of the power efficiency against the V-number is 
dependent on how V is being changed or, in other words, the power 


(Pcore/Prad)*10**5 
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Figure IV-6a. Power efficiency of a bulk distribution of sources versus the V- 
number. This graph was obtained by varying the wavelength, X, from 
0.3Aan to 2.0/zm. The core radius, the cladding radius and the indices of 
refraction of the core and cladding are held fixed at a=5.0/xm, b=25.0/im, 
n, =1.6 and n Had =1.55. 



(Pcore/Prad)*10 
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Figure IV-7a. Power efficiency of a bulk distribution of sources versus the V- 
number. This graph was obtained by varying the wavelength, X, from 
0.45^m to 2.0nm. The core radius, the cladding radius and the indices of 
refraction of the core and cladding are held fixed at a= 5.0/xm, fc=25.0^m, 
m =3.0 and 1.1. 
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Figure IV-7b. Same as Figure IV-7a, using the wavelength, X, instead of the 
V-number. The core radius, the cladding radius and the indices of 
refraction of the core and cladding are held fixed at a=5. 0/im, b= 25.0/tm, 
« core =3.0and« clad =l.l. 
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efficiency does not always increase with the V-number. Later results of this 
work have been consistent with the conclusion that the power efficiency does 
increase with the difference in indices of refraction, n core -n clad a conclusion 
which can explain Figures (IV-1) thru (IV-3). 

4) In the range of small V-number the power efficiency tends to increase 
slightly (Figure (IV-6a)). Such behavior can also be observed in Figure (IV-6b) 
in the high wavelength' region. Maybe, at small V-numbers, the rate of 
increase in the number of modes tends to catch up with the decrease in the 
wavelength. The final result is an increase of the power efficiency at low V- 
number and a gradual decrease for higher V-numbers. We did not observe 
the same behavior in Figures (IV-7) most likely because we could not reach 
the low V-number region (the higher the differences in the indices of 
refraction the larger the V-number). 

ii) Thin Film Distribution 

Figures (IV-8) and (IV-9) reproduce the results of the thin film for a variable 
wavelength. The values used for the parameters of Figures (IV-8) were the 
same as those in Figures (IV-6). We have used the indices of refraction 
n COT =2.5 and n ciad =1 - 4 in Fi 8 ures ^ IV ' 9) ' AU other P arameters unchanged 
with respect to the previous figures. Notice that the P eff remains almost 
constant over the whole spectrum of values. Probably, the tunneling effect 
that occurs at a high wavelength is no longer important for the thin film 
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because of the location and/or reduced number of layers of sources in the 
cladding. 

IV-6 Power Efficiency Versus Core Radius 

Figures (IV-10) and (IV-11) display the behavior of the P eff with the core 
radius for bulk and thin film distributions respectively. For both cases we 
have used the values H cor =1.5 and A=0.6//m. For n clad we used 1.4 and 1.0 
for bulk and thin film respectively. Notice that for the bulk distribution we 
have obtained a general decrease in P eff . The peaks observed are the ones 
corresponding to the cut-off frequencies of the modes. Again, P eff decreases 
with the V- number which varies from 5.64 to 56.4. It should be remembered 
that there is a linear relationship between the core radius and the V-number. 

For the thin film case, we have an increase in the power efficiency, somehow 
abrupt, at low V-number, and a slight decrease at high V-number. It is 
interesting to notice that, for some reason, the decrease of the P eff with the 
fiber radius is more pronounced in the bulk case than in the thin film case. In 
both cases this decrease may be related to the fact that a big core radius inhibits 
the tunneling from the cladding. 

IV-7 Power Efficiency Versus « cor ;« c ] a d (V Constant) 

Figures (IV-12) and (IV-13) illustrate how the power efficiency behaves with 
the difference « C0re -« C | a< j and constant V for bulk and thin film distribution 
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Figure IV-8a. Power efficiency of a thin-film distribution of sources versus the 
V-number. This graph was obtained by varying the wavelength. A, from 
0.3fim to 2.0fj.m. The core radius and the indices of refraction of the core 
and cladding are held fixed at a= 5.0/im, n core =T6 and n c i a d =1 - 55 - 
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Figure IV-8b. Same as Figure IV-8a (thin-film) using the wavelength, X, 
instead of the V-number. The core radius and the indices of refraction of 
the core and cladding are held fixed at a=5.0fim, « core = 1.6 and « clad =1.55. 
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Figure IV-9a. Power efficiency of a thin-film distribution of sources versus the 
V -number. This graph was obtained by varying the wavelength. A, from 
0.3^/m to 2.0/im. The core radius and the indices of refraction of the core 
and cladding are held fixed at a=5.0/im, n cor =2.5 and ^ c i a d =1 - 4 - 
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Figure IV-9b. Same as Figure IV-9a (thin-film) using the wavelength, X, 
instead of the V-number. The core radius and the indices of refraction of 
the core and cladding are held fixed at a=5.0^m, n core = 2.5 and n clad =1.4> 
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Figure IV-10. Power efficiency of a bulk distribution of sources versus the 
radius a. The cladding radius is given by 5.0a. The V-number varied 
from 5.64 to 56.4. The wavelength, A, and the indices of refraction of the 
core and cladding are held fixed at A=0.6/im, n core =1 - 5 and n dad =1 - 4 
respectively. 
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Figure IV-11. Power efficiency of a thin-film distribution of sources versus the 
radius a. The V-number varied from 11.7 to 186.74. The wavelength. A, 
and the indices of refraction of the core and cladding are held fixed at 
A=0.6/im, n core =1.5 and tt clad =1.0 respectively. 
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respectively. The results were obtained for four different V-numbers and 
were plotted in a log-log scale. Using a linear equation we have found a 
correlation coefficient of one for the first three V-numbers. For the forth V- 
number, V=62.83, the correlation coefficient was 0.998 for a linear equation 
and 1.0 for a quadratic one. Apparently, the higher the V-number the more 
the graph deviates from a linear equation in a log-log scale. These show that 
the greater the difference between the indices of refraction, the higher the 
power efficiency. As discussed before, the apparent high value of the power 
efficiency at high V-number is due not to the increase in this value but to 
bigger differences between the indices of refraction. The upper portion of each 
curve is the highest power efficiency that can be reached at these particular 
values. That is because at those points, the index of refraction of the cladding 
is equal to one (the lowest possible index of refraction and the highest possible 
difference between the indices at V constant). As both indices of refraction 
increase, the difference between them and the power efficiency decrease. The 
indices of refraction obey the inequality 3.5>n cor ^>« clad >10. In this way, the 
lower extremes of each curve corresponds to w cor e = 3.5. It may be argued that 
this value of the index of refraction is not realistic. However, as a theoretical 
work, extreme cases like this one, are of interest because they help 
determining the general behavior of this system. Notice that the V-numbers 
used have a very big difference between each other. Have we had used values 
which were closer to each other we could have observed lines that obey the 
classical increase of the power efficiency at the cut-off value. 
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Figure IV-12. Power efficiency of a bulk distribution of sources versus the 
difference ” core -« c i ad - The indices of refraction are restricted to the 
interval 3.5>n cor ^>n c | ad >1.0 under the condition of constant V-number. 
The smaller the difference, the greater the indices of refraction. The core 
radius, the cladding radius and the wavelength. A, are held fixed at 
a=6.0/im, fc= 30.0/im and A=0.6^m, respectively. 
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Figure IV-13. Power efficiency of a thin-film distribution of sources versus the 
difference n core -« clad . The indices of refraction are restricted to the 
interval 3.5>n cor ^>« clad >l-0 under the condition of constant V-number. 
The smaller the difference, the greater the indices of refraction. The core 
radius and the wavelength. A, are held fixed at a=6.0/xm and A=0.6Aim 
respectively. 
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IV-8 Bare Core Fiber with a Thin Film Coating 

We conclude our data analysis with an additional graph that simulates a very 
useful fiber sensor: the bare fiber core coated with a thin film of fluorescent 
material. Such fiber can be modeled as a cylindrical rod whose cladding is the 
air itself (rc c i ad =1.0). The coating is assumed to be thin enough so the effect of 
its index of refraction can be ignored. In order to know what is the behavior 
of the power efficiency of this fiber for a variable « cort? we have plotted 
Figure (IV-14). As it was already expected, the power efficiency increases with 
the "core ( the higher the difference « core -n clad , the higher the power 
efficiency). The increase in the power efficiency is almost linear but can not 
keep on growing indefinitely. Marcuse has found a power efficiency as high as 
60% [Marcuse, 1988]. Our graph shows that it could go even higher, 85%. If we 
could use larger and larger core indices of refraction, we should expect higher 
efficiencies. However, we should point out that the results discussed refer 
only to the forward propagating modes of the fiber. Consequently, if we take 
into account both forward and backward propagating modes, the final result 
would exceed the 100% limit of the power efficiency! Later on, in a personal 
communication to Marcuse, we have suggested that the total power radiated, 
P rad' ma y had been understimated by a factor of two. However, in a 
subsequent communication, Marcuse has stated that although the total power 
radiated may had been understimated, the correct expression for P rad may 
include the sum of both forward and backward propagating bound modes 
yielding the following result for the power efficiency 
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pi Eure jv-14. Power efficiency of a Ihin-film distribution of sources versus the 
index of refraction of the core n co , e The core radius, the wavelength and 

the index of refraction of the cladding are held fixed at a=5.0pm, 
A=0.6/im and M cl , d =1.0 respectively. The V-number varied from 29.73 to 
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(IV-6) 


as opposed to 


Peff = • (III-34a) 

P rad 

So, i/ Equation (IV-6) is correct, the problem of the P eff greater than 100% 
has been solved. Equation (III-34a) would still be a good approximation 
whenever P rad >»P core or P eff is small enough (lesser than 0.1) and most of 
the plots presented, particularly the ones of the bulk distribution, would be 
valid. 

Marcuse has cautioned (in a private communication) that the expression for 
Prad ^ still approximate. More specifically, back in the calculation of P rad , 
Section (III-4), a plane wave was used to describe both electric and magnetic 
fields of the radiation modes [Marcuse, 1975]. Accordingly, the expression of 
the total power radiated, P' rad , has to involve the actual fields of the 
radiation modes. Although these expressions were presented by Snyder et. al. 
[1983], Chapter 25, Snyder [1971] and Sammut [1982], they haven't been used 
yet in the corresponding integral. However, since at infinity the radiation 
fields are well described by a plane wave, it is our belief that the plane wave 
approximation is a good one. 


V. CONCLUSIONS AND FURTHER WORK 


Previous work on the light injection efficiency into the core of a fiber from 
sources in the cladding have made use of the fields of a weakly guiding fiber 
[Marcuse, 1988]. This approximation simplifies the analysis of the power 
efficiency by introducing universal values for the eigenvalues of different 
fibers with same P-number, but can not predict accurately the behavior of 
the injected light into a fiber with arbitrary differences in indices of refraction. 
Although Marcuse has analyzed the behavior of his weakly guiding fiber as a 
function of the P-number, he did not do so for the' other parameters. We 
have extended Marcuse's work by using the exact field solution in the 
expressions of the power efficiency, P eff , and analyzed its behavior as a 
function of the remaining parameters. Although more complicated and 
harder to interpret, our formulas allow us to analyze the power injection 
efficiency of fibers with arbitrary differences in the indices of refraction. The 
results obtained are relevant for the design of more efficient optical fiber 
distributed sensors and are summarized below. 

We have confirmed the weakly guiding results of Marcuse [1988]. However, 
we found that the P eff does not always increase with the P-number but 
with the difference in the indices of refraction, n core^ciad' This conclusion 
could not be easily reached by Marcuse because of his approximation. 
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For a fixed ratio a / X, outer radius of the distribution and indices of 
refraction, the P eff is independent of the core radius, a, and the wavelength, 
X. This suggests that a / X is an independent variable. 

We have also found that, in general, the P e (( of the bulk distribution 
increases with the wavelength, X. However, for the thin film, the P eff 
remains almost constant. The first result implies a decrease of the power 
efficiency with the U-number which is contrary to the general belief that the 
power efficiency would increase with the V'-number [Marcuse, 1988; 
Lieberman et. al., 1990]. 

The behavior of the P eff with the fiber radius is also different for the bulk and 
thin film distributions. For the bulk case the power efficiency decreases for 
values of the core radius. However, for the thin film it increases at 
small values of core radius and decreases slightly for bigger values. Again, we 
have obtained a decrease of the power efficiency with the Wnumber. 

A final result states that the power efficiency increases with the difference 
between the indices of refraction H core -n clad for both bulk and thin film 
distributions. Consequently, in order to build a fiber sensor which has the best 
performance, we should use a source coating that emits light of high 
wavelength, use the smallest possible core for fibers with a bulk distribution 
and a slightly larger one for a thin film fiber. A fiber with a large refractive 
indices differential is optimum for both distributions. Nothing can be 
concluded for the dependence of the power efficiency with the cladding radius 
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b because of the infinite cladding approximation used at the very beginning. 

It is hard to compare a thin film fiber coating with the bulk one. Although, 
the efficiency of the thin film is at least three orders of magnitude higher, 
what really matters is the total signal collected at the end of the fiber [Marcuse, 
1988]. Experimental results by Glass et. al. [1987] on the signal response of both 
distributions have determined that the bulk distribution signal is higher. 
However we must have in mind that, for both distributions, most of the 
evanescent interaction occurs at the core/cladding boundary, for this reason, 
it is desirable to accumulate as many sources as possible in this region. 

Due to the infinite cladding approximation, the fields used to compute the 
power efficiency are independent of the cladding radius, b. The later 
dependency in b in the expression of the power efficiency of the bulk 
distribution arises because of the finite distribution of the cladding sources, in 
other words, the parameter b that appears in these equations is related to 
r out and not to the claddin 8 radius itself! Additional work is required to 

determine a more accurate behavior of the P eff with b. 

The plane wave approximation used to calculate the total power radiated 
could be tested if we could integrate Equation (III-27) using the actual 
radiation modes of the fiber. Apparently, it can be done. However, if that is 
not the case, we still could compare both results (plane wave and the 
expression for the radiation modes) using the special case of a current source 
j polarized in the z-direction. Choosing this polarization could simplify the 
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integration of Equation (III-27). 

Further experimental work should be undertaken to verify our model. This 
seems to be the most logical follow up to this work. Such an experiment was 
conducted earlier by Glass et. al. [1987], Love et. al. [1988] and Lieberman et. al. 
[1990]. The first two authors used a ray theory approach (geometric optics) to 
compare their experimental results. They were mainly interested in the total 
power injected into the fiber. However, Lieberman et. al. [1990], have 
determined the power efficiency of a fluorescent cladding fiber for a single 
experimental data point. They have found an agreement within an order of 
magnitude; a result that lends support to Marcuse's theory and the current 
work as well. 

In summary, the richness of results provided by this rather simple fiber, the 
circular step index profile, suggests that more complicated fiber geometries 
and index profiles may have many other characteristics well worth exploring. 
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